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Abstract 

We investigate the effects of a time-correlated noise on an extended chaotic system. The chosen 
model is the Lorenz'96, a kind of "toy" model used for climate studies. Through the analysis of the 
system's time evolution and its time and space correlations, we have obtained numerical evidence 
for two stochastic resonance-like behavior. Such behavior is seen when both, the usual and a 
generalized signal-to-noise ratio function are depicted as a function of the external noise intensity 
or the system size. The underlying mechanism seems to be associated to a noise-induced chaos 
reduction. The possible relevance of these and other findings for an optimal climate prediction are 
discussed. 
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I. INTRODUCTION 



The last decades have witnessed a growing interest in the study of the effect of noise on 
dynamical systems. It was proved that, under some conditions, when a nonlinear dynamical 
system is subject to noise, new phenomena can arise, phenomena that only occur under the 
effect of such noise. All these phenomena are lump together under the name noise-induced 
phenomena. A few examples are: stochastic resonance in zero-dimensional and extended 
...ten,. II a a. noi.e-i„duced transition. W. noi.e-induced phase transition. Q £ noi.e- 

nn, nn n, 

induced transport [7|, |8[, noise-sustained patterns [9|, llO[, noise-induced limit cycle 

Clearly, some of the above indicated noise-induced phenomena occur in spatially extended 



systems, where another phenomena of great relevance exists: spatio-temporal chaos 
However, studies on the effect of noise on spatially extended chaotic systems are scarce 
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There are studies on chaotic systems where the pseudo-random behavior of the system is 



the trigger o: 
for instance 



lenomena usually associated with the effect of a real stochastic process (see 
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151]). Hence, we can refer to the presence of a deterministic noise, that is a 



pseudo-noisy behavior associated to the chaotic character of the system. 

Among others, one of the most relevant and largely analyzed application of studies of 
chaos in extended systems corresponds to climate prediction. These kind of problems have 
been described by Lorenz |17j as falling into two categories. On one hand those which 
depend on the initial conditions, while on the other are those depending on the boundaries. 
However, both kind of prediction problems are affected by errors in the model equations 
used to approximately described the behavior of real systems. 

In a recent work [16| , and with the aim to improve the limited weather predictability that 
results from a combination of initial conditions uncertainty and model error, it was presented 
a study of the effect of a stochastic parametrization in the Lorenz'96 model 
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19|. It is 



well known that much of the current error in weather predictability derives from the practice 
of representing the effects of process occurring at unresolved scales by using simple forms of 
deterministic parametrization, attempting to summarize the effects of small-scale processes 
in terms of larger-scale, resolved, prognostic variables. 

In this work, and with a similar objective as in a previously indicated paper [loj], we 
investigate the effect of a time-correlated noise on an extended chaotic system, analyzing the 
competence between the above indicated deterministic noise and a real stochastic process. 
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In order to perform such a study we have chosen the Lorenz'96 model 17[_^ In spite of the 
fact that it is a kind of toy-model, at variance with the cases studied in [13], it has a clear 
contact with real systems as it is of interest for the analysis of climate behavior and weather 



prediction 
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191]. In fact, this model has been heuristically formulated as the simplest 



way to take into account certain properties of global atmospheric models. To reach our 
objective, we have assumed that the only model parameter is time dependent and composed 
of two parts, a constant deterministic contribution plus a stochastic one. 

Through the analysis of the system's temporal evolution and its time and space cor- 
relations, we have obtained numerical evidence for two stochastic resonance-like (SR) |l| 
behaviors. Such behaviors are seen when both, the usual signal-to-noise ratio (SNR) and 
a generalized function SNRgiob (that we call global SNR), are depicted as function of the 
external noise intensity or the system size. In accord with what was shown in previous works 
|l|, a SR phenomenon can occur in systems without external periodic forcing, but having an 
internal typical frequency. Hence, it seems reasonable to assume that the present resonances 
typically occur at frequencies corresponding to a system's internal quasi-periodic behavior, 
as well as at an optimal system's size. Finally, we discuss the possible relevance of these 
findings for climate prediction. 



II. THE MODEL AND RESPONSE MEASURES 



The Model Lorenz'96 



The equations corresponding to the Lorenz'96 model 17|, ll8| are 



ij(t) = —Xj^i{Xj^2 — Xj+i) — Xj + F. (1) 

J = 1,2,3,. ..,iV, 

where Xj{t) indicates the time derivative of Xj{t). In order to simulate a scalar meteorological 
quantity extended around a latitude circle, we consider periodic boundary conditions: Xq = 

XNi X^i = XAr_i. ^ ^ 

As indicated before, the Lorenz'96 model ^/3] has been heuristically formulated as the 
simplest way to take into account certain properties of global atmospheric models. The terms 
included in the equation intend to simulate advection, dissipation and forcing respectively. 
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In contrast with other toy models used in the analysis of extended chaotic systems and based 
on coupled map lattices, the Lorenz'96 system exhibits extended chaos {F > 9/8), with a 
spatial structure in the form of moving waves The length of these waves is close to 5 
spatial units. It is worth noting that the system has scaled variables with unit coefficients, 
hence the time unit is the dissipative decay time. In these units the group velocity of the 
waves is close to Vgr = 1.20 implying a eastward propagation. If, as done by Lorenz, we 
associate the time unit to 5 days and the system size of 40 to the length of a latitude circle, 
we have a highly illustrative representation of a global model. If in addition we adjust the 
value of the parameter F to give a reasonable signal to noise ratio (Lorenz considered F = 8) 
the model could be most adequate to perform basic studies of predictability. Hence, within 
this framework, the signal analyzed in this paper would correspond to the passing of waves 
in a generic observational site, in what is a simple mimic of forecasting at an intermediate 
time range. 



B. Stochastic contribution 



As indicated before, here we assume that the model parameter F becomes time dependent, 
and has two contributions, a constant and a random one 



F,{t) = Fmed + '^jit), 



(2) 



with ^E'j(t) a dichotomic process. That is, "^jit) adopts the values ±A, with a transition 
rate 7: each state (±A) changes according to the waitin g ti me distribution ipiit) ~ c 
The noise intensity for this process is defined through |20|, |21| ^ 
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C. System Response 

As a measure of the SR system's response we have used the signal-to-noise ratio (SNR) 
To obtain the SNR we need to previously evaluate S{uj), the power spectral density 



(psd), defined as the Fourier transform of the correlation function 20|, |2l| 
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S{u;)= / e^-^(x,(0)a;,(r))rfr, (3) 
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where ( ) indicates the average over realizations. As we have periodic boundary conditions 
simulating a closed system, (xj(O)xj(r)) has a homogeneous spatial behavior. Hence, it is 



enough to analyze the response in a single site. 

We consider two forms of SNR. In one hand the usual SNR measure at the resonant 
frequency Uo (that is, in fact, at the frequency associated to the highest peak in S{uj)) is 

SNR = (4) 

where 2a is a very small range around the resonant frequency Uo, and Sbacki^^) corresponds 
to the background psd. On the other hand we consider a global form of the SNR (SNRgiob) 
defined through 

J^2l' dcuSitu) 
SNRgiob = f-urZTi e r^' 

J _ dujbback{UJ) 

where uOmin and uJmax define the frequency range where ^(ci;) has a reach peak structure 
(with several resonant frequencies). 



III. RESULTS 



We have analyzed the typical behavior of trajectories as Xi(t) — Xmed-T, where Xmed-r is 
the time average. It is worth commenting that when the Lorenz96 system evolves without 
external noise (that is Fj{t) is constant), the time evolution shows a random-like behavior, 
with the main feature that the amplitude of the oscillator is constant over all the time. 
However, when the system is subject to a random force as described in Eq. (j2j), the temporal 
response decays, due to the fact that the interaction between the intrinsic evolution and the 
external noise produces a dissipative contribution on the system. Hence the system's time 
evolution consists of a transitory regime and a stationary one. This was analyzed through the 
behavior of the "decay" of {xi{t) — Xmed-r)- We assumed that this decay can be adjusted 
by an exponential law. The decay parameter (A) only depends on Fmed and it does not 
depend neither on the system size nor on the noise intensity. This analysis is relevant when 
studying the effects of noise on the stationary regime. From those results it was possible to 
anticipate, and approximately identify, the existence of two regimes: a weak or undeveloped 
chaos for Fmed < 6.0, and strong or completely developed chaos for Fmed > 6.0. 

The typical numbers we have used in our simulations are: averages over 10^ histories, 
and ~ 10^ simulation time steps (within the stationary regime, see later). 

We have evaluated the psd S{u!) in a standard way. Figure [T]-a shows the typical form of 
the psd S{uj) for a couple of values of F^ed {Fmed = 4.5, 7.8) and for a noise intensity A = 0.1 
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FIG. 1: (a) S{u)) for a couple of values of Fmed {Fmed = 4.5,7.8) and for a small noise intensity 
= 5 X 10"^, A = 0.1 and N = 128). (b) Spatial spectrum for N = 256, A = 0.001 and 7 = 10, 
for several values of Fmed- continuous line -Fmed = 5, dashed F^ed = 6, dotted F^ed = 7 and 
dash-dotted Fmed = 8 • 

= 5 X 10""^). The figure shows a rich peak structure within the interval 0.22 < u; < 1.3. 
It is worth to comment that the frequencies associated to the different peaks seems to 
correspond to the harmonics of the main (or first) peak frequency. In Fig. [T]-b, we show the 
form of S{k), the associated spatial spectrum. Here we depict the spectrum for fixed values 
of the system's size (A^ = 256), and noise intensity (A = 0.001 and 7 = 10), and different 
values of Fmed- The independence of the position of the peak (indicating a single spatial 
structure of wavelength k/2'K = 0.2) is apparent. However there is a strong dependence on 
the peak intensity when varying Fmed, from a net peak for underdeveloped chaos {Fmed = 5) 
to a reduced peak for well developed chaos {Fmed = §)• It is worth here remarking that there 
is no dependence (or eventually a very weak one) of this behavior on the noise intensity. 

Figure [2] shows the dependence of SNRgiot -for a space-temporal noise- on A for fixed 
values of and a couple of values of Fmed- It shows a peak for A ~ 6 — 7 x 10"^, that 
corresponds to the fingerprint of the more usual form of SR. The insert shows the same case 
but for SNR. 

The analysis of the dependence of SNRgiob on Fmed have also shown the existence of 
the previously indicated two regimes: a weak or undeveloped chaos for Fmed < 6.0, and a 
strong or completely developed chaos for Fmed > 6.0. Those regimes are characterized by 
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FIG. 2: SNRgiob vs A, for fixed values of = 64 and 7 = 10; and a couple of values of Fmed- 
white circles Fmed = 5, black squares Fmed = 6. The insert shows the behavior of the usual SNR 
for the same cases. 

the existence of well defined peaks in the psd, in the former case, and a less defined peak 
structure in the latter case, as seen in Fig. [1^ . 

It is worth to detach the strong similarities in the behavior of SNR and SNRgiob -which 
becomes apparent when comparing the main Fig. [2] with its insert- indicating that the 
second one is an adequate and more versatile measure to characterize the system's response. 
Hence, due to the clearness in the determination of SNRgiob (compared with the difficulties 
for a correct determination of SNR for large values of Fmed) in what follows we adopt it for 
the system's analysis. 

In Fig. [3] we depict the dependence of SNRgiob on A^, for the case of space-temporal 
noise, for fixed values of A = 0.1 and 7 = 5 = 0.001) and a couple of values of Fmed- 
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FIG. 3: SNRgiob vs A^, for fixed values of A = 0.1 and 7 = 5 = 0.001), and a couple of values of 
Fmed'- white circles Fmed = 5, black squares -Fmed = 6. The insert shows the behavior of the usual 
SNR for the same cases. 

The existence of the peak at ~ 60 for Fmed = 5.0 is apparent. In addition, we observe an 
increase of SNRgiob for large values of A^. However, for Fmed = 6.0, the peak has disappeared, 
as well as the increase with larger values of A^. The presence of the peak at A^ ~ 60 indicates 
a kind of system-size stochastic-resonance (SSSR) [23]. The insert shows the same case but 
for SNR. Again, as indicated above, the nice agreement between the behavior of SNRgiob 
and that of SNR. 

The figures clearly show that the system's response (SNR) is stronger when the system is 
in the underdeveloped chaos range than when it is in the highly-developed chaos one. Our 
results also show that the main resonant frequency does not depend on the noise intensity, 
system size, or correlation rate. 

We want to close this section commenting that the SR phenomena found here looks 
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FIG. 4: S{uj) with and without external noise, for F = 5 and 7 = 10. Continuous hne A = 0, 
dot-hne A = 0.1. 



similar to the so called internal signal SR 22|]. In previous studies it was shown that in 



some systems having an internal typical frequency, SR can occur not only at the frequency 
of an external driving signal, but at the frequency corresponding to the internal periodic 
behavior 22| . Regarding the present mechanism of SR, what we can indeed comment is that 
the increase in the SNR is related not to a reinforcement of the peak high respect to the 
noisy background at a given frequency, but with a reduction of the pseudo (or deterministic) 
noisy background when turning on the real noise. That is, the interplay between "real" 
noise and "deterministic" noise conforms a kind of noise-induced chaos reduction. Figure H] 
shows, for fixed values of F and 7, the behavior of S{uj) in both cases: with (A 7^ 0) and 
without noise (A = 0). The above indicated reduction trend, as the real noise is turned on, 
is apparent. However, the present mechanism is not completely clear so far, and requires 
further studies. 



IV. CONCLUSIONS 



We have investigated the effect of a time-correlated noise on an extended chaotic system, 
analyzing the competence between the indicated deterministic or pseudo-noise and the real 



random process. For our study we have chosen the Lorenz'96 model 



17| that, in spite of 
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;he fact that it is a kind of toy model, is of interest for the analysis of climate behavior 
161 . 119|. It worth remarking that it accounts in a simple way of the spatial structure of 
geostrophic waves and the dynamics of tropical winds. The time series obtained at a generic 
site Xi{t) mimics the passing of such waves, which is in fact a typical forecast event. We 
have assumed that the unique model parameter F, is time dependent and composed of two 
parts, a constant deterministic, and a stochastic contribution in a spatial-temporal form. 

We have done a thorough analysis of the system's temporal evolution and its time and 
space correlations. From our results it is clear that, using two complementary SNR measures, 
a usual and a global one, we have obtained numerical evidence for two SR-like behaviors. 
In one hand a "normal" SR phenomenon occur at frequencies that seem to correspond 
to a system's quasi-periodic behavior. On the other hand, we have found a SSSR-like 
behavior, indicating that there is an optimal system size for the analysis of the spatial 
system's response. As indicated before, the effect of noise is stronger when the chaos is 
underdeveloped. 

We argue that these findings are of interest for an optimal climate prediction. It is clear 
that the inclusion of the effect of an external noise, that is a stochastic parametrization 
of unknown external influences, could strongly affect the deterministic system response, 
particularly through the possibility of an enhanced system's response in the form of resonant- 
like behavior. It is worth here remarking the excellent agreement between the resonant 
frequencies and wave length found here, and the estimates of Lorenz 18|. 

The effect of noise is weak respect to changes in the spatial structure, with the main fre- 
quencies remaining unaltered, but it is strong concerning the strength of the "self-generated" 
deterministic noise. In fact, in such a system and at the resonant frequencies, forecasting 
could be improved by the external noise due to the effect of suppression of the self-generated 
chaotic noise. The detailed analysis of such an aspect will be the subject of a forthcoming 
study |24|. 
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